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Abstract 

We consider a non-resonant system of finitely many bilinear Schrodinger equations with discrete spectrum 
driven by the same scalar control. We prove that this system can approximately track any given system of trajectories 
of density matrices, up to the phase of the coordinates. The result is valid both for bounded and unbounded 
Schrodinger operators. The method used relies on finite-dimensional control techniques applied to Lie groups. We 
provide also an example showing that no approximate tracking of both modulus and phase is possible 
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1 Introduction 

a 

' — ^1.1 Physical Context 

^ The Schrodinger equation describes the evolution of the probability distribution of the position of a particle in the 
00 space. The evolution of the Schrodinger equation can be modified by acting on the electric field, e.g., through the 
action of a laser 

r — . 

^ We will be interested in this paper in non-relativistic and non-stochastic Schrodinger equations on a domain (i.e., 
an open connected subset) D. of R'' that is either bounded or equal to the whole R'^ {d G N). To each equation, we 
associate a Schrodinger operator defined as 

0\ . 

\if{x)) ^ {x^ —A\ir{x)+V{x)\ir{x)), x^Q., 

IT/i where Y denotes the wave function and the real-valued function V is called the potential of the Schrodinger operator. 

The wave function verifies V^^ = 1 ■ We assume moreover that V is extended as +oo on R^ \ Q., so that, in the case 
. 5^ Q. bounded, \\f satisfies the boundary condition = 0. The controlled Schrodinger equation with one scalar control 

is the evolution equation 

= -A\l/{x,t)+V{x)\lf{x,t) + u{t)W{x)\lf{x,t), (1) 

where the real-valued function W is the controlled potential. 

The control function m : [0, T] ^ R is chosen among a set of admissible control functions in order to steer the 
quantum particle from its initial state to a prescribed target. A classical result of Ball, Marsden and Slemrod asserts 
that in general, exact controllability is hopeless (see lH and ID for precise statements and a proof). 

The approximate controllability of one particular system of the type ([Hi has already be proved by Beauchard using 
Coron's return method (see [5 1 and references therein). Approximate controllability results for general systems under 
generic hypotheses were proved later with completely different methods in |i6j|. 
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An interesting question is the simultaneous control of several Schrodinger equations of type (O, describing the 
evolution of several particles, driven by only one control. The most elementary case is the one where identical particles 
are submitted to the same control field (that is, Schrodinger operator and controlled potential are the same for every 
particle), but the initial positions and the targets differs from one particle to another. This problem is often referred as 
control of density matrices by physicists. This case has already been successfully addressed in [61. 

In this work, we study the approximate controllability of a system of density matrices, covering among others the 
following two physicals situations: 

In the first situation, we consider a set of particles of different natures, that is, their Schrodinger operators are all 
different. Initial states and targets can be arbitrary. The control function u is the same for all particles. 

In the second situation, we consider a set of identical particles, with the same Schrodinger operator, but with 
different controlled potentials. This second situation is of particular interest since it is usually easier to select a set of 
identical particles than to ensure that they are all submitted to the same excitation. The control function u is the same 
for all particles. 

Our main result is that, under non-resonnance and connectedness conditions, in the both cases above, whatever 
the sources and the target are, it is possible to chose an arbitrary small piecewise constant control function u in such a 
way that every particle approximately reaches its goal (with an arbitrary precision), and moreover that given any path 
joining the sources and the target, it can be followed (up to the phase) by the particles with an arbitrary precision. Our 
approach extends to the control of the density matrices of the considered systems. 

1.2 Mathematical framework 

We give below the abstract mathematical framework which will be used to formulate and prove the controllability 
results later applied to the Schrodinger equation. The fact that the Schrodinger equation fits the abstract framework 
has already been discussed in iH Section 3]. 

Let U he a subset of R. Let p be a nonzero positive integer. For I < i < p, let //, be an Hilbert space, be an 
integer. A,- : D{Ai) C Hi — > Hi be a densely defined (not necessarily bounded) essentially self-adjoint operator and for 
1 < 7 < ni, let Bi j : D{Bij) C Hi —>■ Hi be a densely defined (not necessarily bounded) linear operator. 

We assume that ((A,)i<,<p, (B,j)i<Kpj<y<„., ?7) satisfies the following thi^ee conditions: (HI) A,- and B, y are 
skew-adjoint, (H2) there exists an orthonormal basis (0f of Hi made of eigenvectors of A,-, and all these eigen- 
vectors are associated to simple eigenvalues (H3) (pj^ G D{Bi^j) for 1 < / < 1 < 7 < A crucial consequence 
of these hypotheses is that for every \ <i<p, \ <i<ni and every u ^ U, Ai + uBi j has a self adjoint extension 
on a dense subdomain of //, and generates a group of unitary transformations j) : Hi Hi. In particular, 

^t{Ai+uBij) ^^.-j _ ^, eygj-y u&U,\<i<p,\<j<ni and every t >0, where 5,- is the unit sphere of 

We consider the Yli=\ conservative diagonal single input control systems 



dt 
dt 



Al(V^l,2(0)+w(0^1,2(V^l,2(0) 



< 



{t) 



Al(V^l,„, (0) + (V^l.m (0) 

A2{\lf2,i{t)) + u{t)B2s{\ir2At)) 



(2) 



dt 

dt 



dt 



.np_ 



MVp,n,Xt)) + u{t)Bp^„^X¥pMp (0) 



with initial conditions to be specified later. 

A point 1//° = (Yi itVi2' - ■ ■ ^Vi n,^V2\^ - ■ ■ ^Vp.n,,) of H = nf=i (^;)"' and a piecewise constant function u : 
[0,r] —)■ U, u = Yj^i X[ti,t,+i)^i being given, we say that the solution of ^ with initial condition xj/^ G H and cor- 



responding to the control function u : [0,T] U is the curve 



Y: [0,T] H 

t ^ \lf{t) = {\l/i^i{t),Yl,2{t),---,¥p,n„{t)) 



defined by 



o... oe''^^-+"''^--j\Y?j), (4) 

1 < i < 1 < 7 < by where U <t < L/Li fi and m(t) = uj if Y^jzl ti < r < Y^j^-^ U. Notice that such a Yij{-) 
satisfies, for every « G N and almost every t € [0, T] the differential equation 

A piecewise constant function m : [0, T] ^ R being given, the propagator of the control system in Hi 

=Ai{^^fi,J{t)) + u{t)BiJ{^^fi,j{t)) 

will be denoted by <I>, j . By definition, 

for any t in [0, T] and any \(rfj in 

The propagator of the whole system Q is denoted by <I> = nf=i Il"j'=i ^i j- 

For !</</?, 1 < 7 < «!, and (/:,/) G N^, we define also the numbers ai{k,l) = {Ai^l,^l) and bi j{k,l) = 

A finite sequence {k\,k2, ■ ■ ■ M) of N is said to connect the two levels k and k' for the diagonal conservative 
diagonal single-input control system (Ai,Bij) if ki = k, ki = k' and Y['ij=\ bi,j{kq,kq^i) ^ 0. 

A subset S of is called a connectedness chain of {Ai,Bij) if for arbitrary large N, for A;') in N^, I <k,k' <N, 
there exists a finite sequence {qi,q2, ... ,qi) in [0,N] that connects k and ^' for {Ai,Bij) and such that (^^i^r+i) belongs 
to 5 for every 1 < r < / — 1. 



1.3 Main result 

Theorem 1.1. Assume that the system ([2]) satisfies (i) the concatenation of the spectrum ofAi, A2, Ap is a Q- 
linearly independent family (ii) for any 1 < /q < p, for any I < jo < «,(,, there exists a connectedness chain S of 
(AjQ , Bigjg ) such that for any s in S, for any j in [ 1 , ni^ ] \ { Jo }. | big jg {s) \ ^ \ bi^j {s) \ and ( Hi) U contains a neigborhood 
of zero. 

Let c : [0, r] Y\1^^L{Hi.,HiY' be a continuous curve such that c(0) = Idn and c = IT^i n'yLi curve 
taking value in the set of the Hermitian operators of H such that for every \ <i < p, 1 < 7 < ni, for every time t, 
Cij (t) : Hi — > Hi is unitary in Hi. Let N be an integer Then, for every £ > 0, there exist Tu > T and a piecewise 
constant control u : [0,r„] — > U, such that the corresponding propagator <I> : [0,r„] x H ^ H of system (|2]) satisfies 
i) for every t in [0,7^J, there exists s in [0,T] such that | — |(c,-,y(f)(0/),0^)|| < £ for every k in N, 

^<i<pA<i<ni,\<l<N, and ii) |l<I>ij(r„, 0/) - c,j(r)(0/)|| < efor every I < i < p, I < j < ni, I < I < N. 



1.4 Content of the paper 



To prove Theorem 11.11 we will use a method introduced for the control of the Navier-Stokes equation in IH and 
adapted to the Schrodinger equation in [6|. 

In Section |2l we explain how to choose a Galerkyn approximation of the original infinite dimensional control 
problem ^ in some space SU (k). In Section[3l we use the Lie group structure of SU {k) to compute the dimensions of 
some Lie subalgebras of s\x{k) and to prove that the Galerkyn approximations obtained in Section |2] have some good 
tracking properties. The proof of Theorem 1 1.1 1 and an estimation of the L} norm of the control are given in SectionlH 
where we prove that the original system (O share the ti^acking properties of the Galerkyn approximations established 
in Section[3] A partial counterpart of Theorem I 1 . 1 1 (impossibility of approximate tracking of both the modulus and the 
phase) is given in Section [5] 



2 Choice of Galerkyn approximations 
2.1 Control and time-reparametrization 

We may assume without loss of generality that U has the special form U =]0, 5]. Remark that, if m / 0, = 
g/M((i/«)A,+B,j) Associate with any piecewise constant function u = ^4=1 X[ti.ti+i)'^i in PC{[0,Tu],U) the function v = 
Lf=| Z[t,,t,+i[1/"/ ^ PC{[0,Ty],l/U), with Ty = -ti) and T/ defined by induction by Ti = ti and T/+i = 

T/ + uj{ti^i — ti). Up to the time and control reparametrization given above, it is enough to prove Theorem 1 1.1 1 for the 
system ((Bi,;),=i..pj=i..„, , (A,-),-=i..p, [^,+°°[): 



dt ' 
dt ' 



v(0Ai(vai,i(0)+6li('/^i,i(0) 

V (OA 1 ( l/Ai ,2 (0 ) + ^1,2 ( V^l ,2 (0 ) 
v{t)Ap{\lfp^n„ (0) + Bp,n,XVp,n,, (0) 



(6) 



where the set of admissible controls is the set PC (R+ , , +00 Q . 

Remark 2.1. A feature of this reparametrization of the control from u : [0, r„] U to v: [0, Ty] 1 /U is that ||m||/^i 



2.2 Galerkyn approximation 

For a fixed piecewise constant control v : R+ -^\/U and a fixed family (v/fy)(=i..p,;=i..«, inH = Y[f=i H"', we consider 
the solution {Yi.j)i=i..p.j=i..ni : t H of the system (O of conservative diagonal single-input control systems with 
initial conditions iw?,j)i<i<p,i<j<nr 

For 1 < / < 1 < 7 < «i, ^ € N , we define the function = {\ifij,<^f) : R ^ C. With our definition of solution, 
x'l j is absolutely continuous and for almost all t in R+ 

= v{t)ai{k,k) + £ bij{k,l)x[j. 

Proposition 2.1. For every \ <i < p and any continuous curve s : [OjTv] Hj taking value in the unit sphere of Hi 
(that is, 11^(011 = ^ for all t in [0,r^]j, define the family f = \{s,<^l)\^, / G N. Then, for any strictly positive £, there 
exists an integer N{£) such that for all t in [0, T^], Y^fl^^ fi{t) > 1 — £. 

Proof. Define gi = fk- For all integer /, gi : [0, Tg] — > R is a continuous function. For any e, the set ^/(e) = € 
[0,Tc]/gi{t) > 1 — £} is an open subset of [0,Tc] (for the topology induced on [0,7^] by the Euclidean distance in R). 
Since for all t, gi{t) is a non-decreasing function of /, it is also clear that for any £, ^/(e) C {e). 



Since, for all t, gi{t) tends to 1 as / tends to infinity, we can say that for all e. 



U/eN^/(e) = [0,r,]. 

Using the compactness of [0, Tc], we can find a finite set Zi , /2, • • • h of integer such that 
Define now A^(£) = sup(/i , Z2, . . . , Ip). For all t in [0, T^], ifif fk{t) > 1 - £. 



□ 



We define nf : Hi C" by 7r,'"(v) = I^'Lj (v, 0f for every v in where Ck is the ^''^ element of the canonical 
basis of C". 

Proposition 2.2. F/jc a reference curve c : [0, T] — > n^Li L{Hi,Hi)"' as in the hypotheses ofTheorem \l.l\ £ > and N a 
positive integer Then, for every {i,j), there exists a continuous curve Mi j : [0, T] —^ SU (m) such that \\n'"{cij{t)^^) — 
Mij{t)Kf(^^\\ < £ for every t in [0,T] and every k in {l..A^}. 



Proof. For I <i < p, \ < j <ni, \ <k <N , apply the Proposition 12. ll to the curve s:t \ 



[t){^l)to get some integer 



m 



For any integer m > sup{m,j i-, I < i < p,l < j < ni, \ < k < N}, for every t in [0, T^], for every I < i < p, 



\<j< ni and every l<k<N, ||nf {{cij{t){\lff)) 



(OW)ll/<£- 



The application Fl^'j : [0,T] (C'")" defined for every t in [0,r] by F!"j{t) = in^jicij(l>^))k=i..N associate with 



every t a family of vectors of {C")^ that is almost orthonormal. Actually, denote ^ : (C")^ (C'")'^ the map that 
associates to any linearly independent family (vi,..,va?) of A'^ vectors of (C'")^ the orthonormal family J^(vi, ..,vn) 
deduced from (vi,..,VAr) by the Schmidt orthonormalization procedure. For any / e {1../?}, j € {!..«/}, t G [0,T], 
F;"j{t) - yfl^jit) tends to G {C"f as m tends to infinity. Noticing that E'/^j : t ^ \\Fl"j{t) - yFI"j{t)\\ is continuous, 
considering the open set J2,„{e) = {? G [0, T] \Efj{t) < e} and applying the same compacity argument as in the proof 
of Proposition 12. 1[ we find an integer m,j- such that — ^fI"-'^ {t)\\ < e for every t G [0, T]. Choosing m = 

suPi<Kp,i<;<„,"J!,; we have ||(7r,'"(c,-j(O0f))/t=i..iv - ^^^'"(011 < ^ for every t in [0,T], every / in {l..p}, every j in 
{I. Mi}. To conclude, we choose Mij{t) in such a way that its first A'^ columns coincide with the first N elements of 



(0- 



□ 



For I < i < p, I < j < ni, define the rxr matrices AJ" = [ai{k,l)]i<i(i<r and B'^ j = [bij{k,l)]\<kj<r- 
The Galerkyn approximation of order m of the system ([6] is 



dt ' 



dt 



v(OA'r(xu(o)+B'r,i(^u(0) 

v(0A'r(xi,2(0)+fiT2(^l,2(0) 



The system ^ defines a control system on the differentiable manifold nf=i IfjLi ^" 
it is possible to lift it to the group of matrices of the resolvent, 



(V) 



Since the system (|7]) is hnear. 



dxi^i 
~dt 
dx\2 
dt 



(0 



dx, 



dt 



v(OA'r(xu(0)+s'{;i(xu(0) 

V(0A™(X1,2(0) + «S-2(X1,2(0) 
v{t)A^{Xp^n,{t))+Bl„^X^p^n,{t)) 



(8) 



which a system in nf=i n"Li with matrix unknowns (xi 1 ,xi,2, -.jXp 

We now proceed to a technical change of variable (variation of the constant) and define for every piecewise 
constant control function v in PC(R, l/U), every positive t and every 1 1 < 7 < jjj- 



Recalling that for all m x m matrices a,b, e "be" = e"'^''"'^b, one checks that (3'(j)i<(<p,i<y<n, verifies 



dy\.2 (^\ 



dt 



dyp. 



np 



dt 



^adiJ^v{s)dsA'l')gm^y^ ^ (f) 
^adU^v{.s)dsA",')gm^^^^^f^ 



(9) 



The system ^ defines a control system on the differentiable manifold n^Li ('«)"'> and for every positive t, every 
l<i<p, every 1 < 7 < and every l<k<m,\ {xij{t),<^f) \ = \ (j,- y(f), 0f ) |. 

Using the canonical injection of nf=i {m)"' in SU {mY^f^y «,), we can write the system Q as 



dy 
dt 



(10) 



where : R — > su (wj^f^i '^i) is ^ that associates to any t the (nil^f^i «,) x (ni^^f^j «,) diagonal block matrix 
constructed from (Lf=i"i) x (L^Li'^O blocks of size mxm, and whose diagonal is (e"'''*'"^"*'Bf"^)i<i<p,i<y<n,- 



3 Tracking properties of the Galerkyn approximations 
3.1 Finite dimensional tracking in semi-simple Lie groups 

First, we have to recall some classical definitions and results for invariant Lie groups control systems. The now 
classical control extensions techniques used here have been introduced by Kupka in the 70's, see |i8.] and references 
therein for details. Most of the material below can be found with proofs and references in 

Let G be a semi-simple compact Lie groups, with Lie algebra g = Tj^G and Lie bracket [,]. The Killing form (see 
l?! Chapter III] for details) is negative definite on the Lie algebra g of G. Its opposite (, ) is a scalar product on g that 
can be extended on the whole tangent bundle TG by the left action of G over itself. The smooth manifold G endowed 
with this bi-invariant scalar product (, ) turns into a Riemannian manifold, whose distance is denoted by do- 

Consider a smooth right invariant control system on G of the form 

U(0) = go ^"^^ 

where ?7 is a subset of R, m : R — > ?7 is a control function to be chosen is one of the following class of regularity 
= { absolutely continuous, measurable bounded, locally integrable, piecewise constant }, / : ?7 ^ g is a smooth 
application, ^0 is a given initial condition and dRab denotes the value of the differential of the right translation by a 
taken at point b.lfG = SU (n), [,] is the standard matrix commutator, the exponential map exp : g — > G is the standard 
matrices exponential, the elements a of G are the matrices with determinant equal to one and such that 'da = Id and 
the elements ^ of g are the zero trace matrices such that 'b + b = 0, and dRab = bxa where x is the standard matrices 
multiplication. 

We define the set y = com {f{u),u € ?7} as the topological closure of the convex hull of all admissible velocities 
at point Id. 

It is obvious that the topological closure of the convex hull of all admissible velocities at point g is dRg{Y). 
We will need the following pretty standard relaxation result: 

Proposition 3.1. Let P be a Lie-subgroup ofG with Lie algebra p.IfY contains some bounded symmetric set S such 
that p C Lie{S), then for any continuous curve c : [0, T] —i- P, for any e >0, for any regularity class k of .J€ , there exist 
Ti, > 0, a control function u : [0,Ti,] U of class k, and an increasing continuous bijection : [0,r„] [0,T] such 
that the trajectory g : [0,Tu] G of {JL) with control u and initial condition c(0) satisfies (i) dG{c{^{t)),g{t)) < sfor 
every t in [0, T,,] (ii) ^{Tu) = T and (Hi) c{T) = g{Tu). 



Proof. Fix e > and a continuous curve c : [0,r] — > P. Up to translation by c(0), one may assume without loss of 
generality that c(0) = Idc- By a classical density argument, one may also assume that for all t, c{t) = exp(v(f)) where 
V : [0,r] ^ p is a piecewise constant function. It is hence enough to study the case for which c{t) = exp(fv) for any t 
with some constant v in p. 

Define by induction 5''"' = S and 5^'+^^ = spanS+ [S,S'] for any integer /. By hypothesis, 5°° = U,£n5'^'^ D p. We 
proceed by induction on / to prove the result if c has the special form c{t) = exp(fv) for any t with some constant v in 

The case where v is in 5 follows from IS Theorem 8.2]. 

Fix / in N and assume that the result is known for any element in S^'\ Choose v in 5^'+^^ and write v as a the limit 
of a (fixed, this is a consequence of the theorem of Caratheodory) linear combination of brackets of elements of S and 

v=lim£A,-[v7,v'7] 

lel 

where / is a finite set. A,- is a (constant) real number, (v'")„ is a sequence of elements of 5, converging to some vj 
and {v"'l)n is a sequence of elements of S^'\ converging to some v". Using once again the symmetry of S and a time 
reparametrization, one may assume that < A,- < 1 for every / and Y,iei = 1- 
Recall now the Baker-Campbell-Hausdorff formula (see Q): for any u, v in q, 

e\p{—tu)exp{—tv)exp{tu)e\p{tv) = exp (t^[u,v] + oo{t^)) , 

that is JG(exp(— fM)exp(— fv)exp(fM)exp(fv),exp (t^[u,v]) = t^a{t), for some function a : R ^ R with limit zero at 
zero. For Tq small enough to be fixed later, define the 4To-periodic piecewise constant function F : R — > g by F{t) = v 
for < f < To, F{t) = u for To < f < 2to, F{t) = — v for 2to < f < 3To and F{t) = —u for Stq < f < 4to, and consider 
the curve g'.t^ exp{F{t)). For any « G N, for any t in [0,?iTo], dG{g{t),Qxp {f-[u,v\) < JJk=i '^o^i'^o) = nT^a{zo). 

Fix T] > and T > 0. Choose Tq small enough such that o;(to) < ^ and n = ^. One gets dG{g{t),exp (f^[M, v]) < rj 
for any t such that <t <T. Apply this last inequality with m in S and v in 5'. The proof of Proposition 13.11 follows 
from [2, Theorem 8.2]. □ 

To obtain trackabillity properties for the system it is enough to check that the finite dimensional systems ^ 
satisfies the conditions on 5 given in Proposition 13 . 1 1 for a suitable p. 

We define the set f = conv({^v(f),v G PC (R, i-^, +°°[) ,t G R+}). In the sequel of this Section, we prove that 
under the hypotheses of Theorem ll.il it is possible to find a set S inY satisfying the hypotheses of Proposition 13. II 



3.2 Some Lie algebraic methods 

Fix an integer m in N, and, for 1 < i < p,l < j < rij,! < k,l < m, denote with Eij^j^j the square matrix of order 

m (r^q^i fiq^ whose entries are all zero but the one with index (m(/ — 1 ) + m(7 — 1 ) + i,m{i — \ ) +m{j — I) + j) which 

is equal to one. (We consider Ejj^kj as a block-matrix. The two first indices stands for the mxm block, the two 
last indices {j,k) stand for the position of the non-zero entry inside the mxm block of index (/, j).) 

Proposition 3.2. Fix r in N. If a and b are two matrices of {su{r)) such that a is diagonal with Q-linearly inde- 
pendent spectrum and b has entries bk^i, for 1 < k,l < m, then all the matrices bk^iEk^i + bi^kEi^k belong to the set 
conv{{Ad^,^^jr^,^)b;v G PC (R, +-[) ,t G R+}). 

Proposition 3.3. Fix r in N. /fa and b are two matrices of (su (r) ) such that a is diagonal with Q-linearly independent 
spectrum and b has entries bj^^, for 1 < j,k < m, then for every 6 in R, all the matrices e'^bj^uEj^k + ^ '^b^jE^j 

belong to the set conv{{e'"^ ^o''"b;v G PC (R, +°°[) ,f G R+}). 

Proof. The proof can be found in 131 Appendix A]. □ 



Applying Proposition l3.2l to the set defined in Section 3. 1, one gets that all the matrices YFi=\ L"Li l^iji^^ i)Eij,kj 
+ bij{l,k)Ei j i ji belong to y. We define 5 as the set of matiices S = {±Lf^i LyLi i>i,j{^<^J)'Eij,kj +b,- j{l,k)Ei j i k, 1 < 
k,l < m}. Proposition 13.31 (applied with 6 = n) proves that S is actually in Y. By definition, S is symmetric and 
bounded. What remains to prove now is that the Lie algebra generated by S is equal to p = Of^i su(m)"' . 

3.3 Reduction to SU {mp 

Proposition 3.4. Choose any 1 < leqio < p, 1 < 70 ^ nj^, 1 < kQ,lo <m such that {ko,lo) is in connectedness chain of 
{Aig,Bigjg) and bi^jikoJo) ^bi^j^ikoJo) for every j jo. Then the matrix bi^j^{kQ,lQ)Ei^,j^^ko,k +bio,h{h,h)Eiojo,k.h 
is in Lie{S). 

Proof. Since the matrices A, have Q linearly independent spectrum, it is enough to apply Proposition l3.2l to see that 
for any I <iQ < p, \ <k,l < m, every matrix bioj{ko,k))'Eigj^koJQ+bio,j{k,ko)Eigj^i^^ko is actually contained in 
y, hence in S. If ni^ = 1, the Proposition I3.4l is proved. 

If > 1, the Proposition|32]is not enough to guaranty that the matrix Z7;ojo('^o,^o)'£'/ojo,<:o,/o +^/o,;o(^o,^o)£';o,7o,fo,<:o 
is actually contained in Y. Nevertheless, one can prove that these matrices are contained in the Lie algebra generated 
by S. 

Indeed, by Proposition [321 a = '£"tibigj{ko,lQ)Eigj^ko,h +l^hjik,ko)Eigj^ko,io belongs to Y. Using Proposition 
|33]with = f , one gets that b = L"=i i^,oj(^o,^o)£,-o.;.*:o,/o - ibiojOo,ko)Eigj,ko,io- Compute [[a,b],b] = -4 l^'j'li 
\bkji'<^o,k)\'^ {bioji^o,k)EiojMoJo+biojiio^'^o)EiojMoJo)' and by induction ad^^ j^^b = (-1)*^2*+^ l^jl^ \bigj{ko,lo)\^'' 
{bioj-{ko,lo)Eigj^ko,k + t>ioj-{lo,ko)Eigj^ko,io) for every ^ G N (see IS for details). A classical Vandermonde argument 
on the linear independence of the vectors {\biQj{ko,lo)\'^)i<j<ni^ gives the result. □ 

The fact that Lie{S) = p follows from Proposition 13 .4l by the hypothesis of connectedness (see Proposition 4.1] 
for a detailed computation). 



4 Infinite dimensional tracking 
4.1 Tracking in the phase variables 

For the proof of Theorem 1 1.1[ we follow the method introduced in 161. From the application c : R ^ L{H,H) and the 
tolerance e given in the hypotheses of Theorem ll.il we use the results presented in Section II.B to find an integer m, 
the finite dimensional control system Q and the trajectory t Yli.jMij{t) to be tracked in SU ('w^f^^ «,). Proposition 
I3.1| gives the existence of some time Ty > and some control function v in PC{ [0,T^],l/U) such that the corresponding 
trajectory (3^1,1 ,..,yp^n^) of ^ tracks the trajectory 1 1-^ YlijMij{t) with an error less than e on each coordinate. 

Since for every I < i < p, I < j < ni, I < k < m, the sequence {bij{k, /));>i is in i^, there exists some Ni in N such 
that YT=N +1 l^i,;(^)OI^ < 'Wr every I <i < p, \ < j <ni,\ <k <m. The next result asserts that any trajectory of 
the system (flOl ) can actually be tracked (up to e), with the A/^i-Galerkyn approximation of system 

Proposition 4.1. There exists a sequence {vk)k in PC (R+ , ] ^ , +00 [) such that for every I < i < p, 1 < 7 < the 

sequence of matrix valued curves t ^ e J > B)^ converges in the integral sense tot ^ 

where 1 1— > G{t) is some continuous curve in U {Ni — m). 

Proof. The proof is a direct application of JH Claim 4.3], dealing with the convergence of the sequence e"^ ' . 

□ 

Proposition 4.2. For k large enough, the control function v = Vk given by Proposition 14. 1 1 satisfies the conclusion (i) 
of Theorem \1.1\ 

Proof. This is a direct application of Claim 4.4]. □ 



UuMijit) 




—m,m 


G{t) , 



4.2 Final phase adjustment 



After time reparametrization, we get a control function u G PC([0, Tu],U) from v. Up to prolongation with the constant 
zero function, the control function u : [0, r„] — > U obtained in Proposition |4.2| can always be assumed to satisfy r„ > T 
(the prolongation obviously still satisfies conclusion (i) of Theorem ll.il ). 

To achieve the proof of Theorem 11. 1[ one has to change u in such a way that it satisfies the conclusion (ii) of 
Theorem ll.il One gets the result with a straightforward application of E Proposition 4.5]. 

4.3 Estimates of the -norm of the control 

Combining the Remark IZT] and the estimates of (3'. Prop 2.7-2.8], one gets an easily computable estimation of the 
L^-norm of the control u. We denote with IJ.ij{t) = \J {M^l the velocity at time t of the 

trajectory to be tracked in SU (m). 

Proposition 4.3. In Theorem \l.l\ one can choose the control u in such a way that 

~ min,-^y^o<i,/<Ar, \bij{k,l)\ 

This estimate is valid for every {bi_j)jj, yet is sometimes trivial or too conservative when some bj j{j,k) is close 
to zero. For these anisotropic situations, when some directions are much easier to follow than others, one can obtain 
sharper estimates using ^ Theorem 2.13], the expressions being slightly more intricate. 



5 To track both the phase and the modulus is impossible 

In this Section, we give a partial counterpart to Theorem I 1.11 Indeed, we exhibit an example for which it is not possible 
to track both the phase and the modulus. The proof can easily be extended to a wide range of systems. 
Consider one single control system in an Hilbert space H 

X = Ax + uBx 

x(0) = (/.I ^^^^ 

where A : // ^ // is a diagonal operator in the Hilbert base {<^i)ieH of ^> with purely imaginary eigenvalues 
and S is a skew adjoint operator whose domain contains <pi for every / in N, satisfying bij = {B^t, ^j) € R for every 
i,j in N. Define as admissible control functions all piecewise constant functions u -.K R+. For Z G N, we note 
xi = {x,^i) the component of the solution of system (ITTI ) and we define aj = S{{xi), bi = 3 (x;). 

Remark 5.1. In the case where B is bounded, it is possible to define solutions of Alii for u in L} (R,R+). The result 
and the proof below are easily extended to integrable controls that are not necessary piecewise constant (in particular, 
to controls that may be not essentially bounded). 

Proposition 5.1. Assume Ai,Z>2,i > 0. Then, for e < for every piecewise constant control function m : R — > 

R+, there exists T > 0, there exists i in N, / > 1 such that |x,(t)| > £. 

In other words, it is not possible to track with an arbitrary precision the constant trajectory xi = \. 

Proof. We proceed by contradiction and assume that there exists some admissible control function m : R ^ R+ such 
that the corresponding trajectory of ([TT]) remains £-close to 0i for every time. From system (fTTI) . we see that 

d /+" \ 

—xi = iX]_xi +u\^Yj i^^i ' ] ' 



that is 



ai = -AiZ7i+m9? (^£(B0i,0y)xyj , (12) 



h = Aiai+M3 (^£(B0i,0y)xyj . (13) 

For any positive t, the integration of (fT3] ) on [0,t] yields bi (t) = Ai J^ai {s)ds+ /q u{s)YX=2^ii^ii^)'^^^ "^hat is 
-e\\B^i\\ u{s)ds < / u{s)Ybiibi{s)ds = bi{t)-Xil < e - X\{1 - e)t 

Jo Jo -^0 

and 

^' Ai(l-e)f 



£||B</>i|| 

Integrating now d2{s) = —^2^2 (■5') + u{s)Y,i-^2b2jai{s) on [0,f] for any t > 0, one finds 

a2{t)>—X2£— [ u{s)ds\\B^\\e + bj I [ u{s)ai{s)ds > {b2i{'i- — e) — s\\B(j)2\\) [ u{s)ds. 
Jo 'Jo ' Jo 

For £ small enough, K = (^2,1(1 — £) — £\\B(j)2\\) > 0, and from ([T4l) . we get a2{t) > Kt for every positive t. Hence, 
a2{t) tends to infinity as t tends to infinity, what is impossible since \a2\ < ||x|| which is constant equal to 1. This gives 
the desired contradiction. □ 
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